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ABSTRACT 

Self-similar and semi-analytical solutions are found for the height- averaged equations govern the 
dynamical behavior of a polytropic, self-gravitating disk under the effects of winds, around the nascent 
object. In order to describe time evolution of the system, we adopt a radius dependent mass loss 
rate, then highlight its importance on both the traditional a and innovative j3 models of viscosity 
prescription. In agreement with some other studies, our solutions represent that Toomre parameter is 
less than one in most regions on the /3-disk which indicates that in such disks gravitational instabilities 
can occur in various distances from the central accretor and so the /3-disk model might provide a good 
explanation of how the planetary systems form. The purpose of the present work is twofold. First, 
examining the structure of disk with wind in comparison to no- wind solution; and second, to see if the 
adopted viscosity prescription affects significantly the dynamical behavior of the disk-wind system. 
We also considered the temperature distribution in our disk by a polytropic condition. The solutions 
imply that, under our boundary conditions, the radial velocity is larger for a-disks and increases as 
wind becomes stronger in both viscosity models. Also, we noticed that the disk thickness increases by 
amplifying the wind or adopting larger values for polytropic exponent 7. It also may globally decrease 
if one prescribe /3-model for the viscosity. Moreover, in both viscosity models, surface density and 
mass accretion rate reduce as wind gets stronger or 7 increases. 

Subject headings: accretion, accretion disks — stars: winds, outflows — gravitation — quasars: general 



1. INTRODUCTION 

Accretion disk, a system where the disk feeds the cen- 
tral object through accretion, under the effect of viscous 
forces, is believed to be present in very different contexts 
and over a wide range of physical scales. In particular, 
such systems are found in active galactic nuclei (AGN), 
stellar black holes (X-ray binaries), young stellar objects 
(YSOs) and quasars (QSOs). Since the angular momen- 
tum is strictly conserved, there needs to be a process that 
transports angular momentum away to prevent it from 
accumulating on the central object. The main accept- 
able possibility for this is the viscosity inside the disk 
(see Frank et al. 2002). But the molecular viscosity is 
inadequate to transport angular momentum in the disk 
and describe luminous accretion disks, so some kind of 
turbulent viscosity is required. 

Many investigators adopt the so-called a- viscous model 
introduced by Shakura (1972) and Shakura & Sanyeav 
(1973) that gives the viscosity (y) at any radius (r) of the 
disk as a product of disk pressure scale height (H), the 
velocity of the sound (c s ), and a parameter a that con- 
tains all unknown physics. Despite a number of success- 
ful applications of the a prescription, this parametriza- 
tion suffers from a number of inconveniences (see Ab- 
bassi et al. 2006, hereafter AGS06). Although the a- 
prescription is based on a kind of turbulent viscosity but 
there is no direct physical evidence for this as the origin 
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of turbulence. 

One of the most significant agents which may influence 
the equilibrium structure and dynamical evolution of al- 
most any kind of disk, is self-gravity. It is not merely a 
matter of forming systems. In the case of accreting super- 
massive black holes, for instance, during most of the evo- 
lution, they are self-gravitating. Historically, accretion 
disk theory has concentrated on the non self-gravitating 
case, the effects of self-gravity having only been discussed 
occasionally (e.g. Paczynski 1978; Kolykhalov & Sun- 
yaev 1979; Lin & Pringle 1987, 1990). For simplicity, 
traditional models of accretion disks assume geometri- 
cally thin configuration and neglect the self-gravity of 
the accreting material, which signifies that only pressure 
and the gravitational force of the central object support 
the vertical structure of the disk. The study of disk self- 
gravity has been the subject of considerable attention 
in recent years. It can be partly due to improved ob- 
servations, which have shown that in several observed 
systems the disk mass can be high enough to have a dy- 
namical role on all scale disks, from AGN to protostars, 
and partly due to the increased computational resources, 
which have allowed a detailed numerical investigation of 
the development of gravitational instabilities in the non- 
linear regime (Lodato 2007 and references therein). 

The structure of self-gravitating disks has been studied 
both through the self-similar solutions assuming steady 
and unsteady state (Mineshige & Umemura 1996; Mi- 
neshige & Umemura 1997, hereafter MU97; Mineshige 
et al. 1997, hereafter MNU97; Tsuribe 1999, hereafter 
TT99; Bertin & Lodato 1999, 2001; Shadmehri & Kha- 



2 



ABBASSI, NOURBAKHSH & SHADMEHRI 



jenabi 2006; AGS06; Shadmehri 2009, hereafter MS09), 
and through direct numerical simulations (Igumenshchev 
& Abramowicz 1999; Stone et al. 1999; Torkelsson et al. 
2000; Gammie 2001; Rice et al. 2003, 2005, 2010; Rice 
& Armitage 2009; Cossins et al. 2010; Meru & Bate 
2011). Of these, a little attention has been paid to the 
polytropic flows (e.g. MNU97; AGS06). 

Some laboratory experiments of Taylor- Couette sys- 
tems (e.g., see, Richard & Zahn 1999; Hure et al. 2001) 
indicate that, although Coriolis force delays the onset 
of turbulence, the flow is ultimately unstable to turbu- 
lence for Reynolds numbers larger than a few thousand. 
In all kind of self-gravitating disks the Reynolds num- 
ber is extremely high, so it was thought that probably 
the hydrodynamical driven turbulent viscosity based on 
critical Reynolds number has a significant role in the re- 
distribution of angular momentum in the self-gravitating 
accretion disks. The resulting hydrodynamically driven 
turbulence would then transport angular momentum effi- 
ciently Duschl et al. (2000) have proposed a generalized 
viscosity prescription based on the hydrodynamically 
driven turbulence at the critical effective Reynolds num- 
ber, /3-prescription, which can be applied for both self- 
gravitating and non self-gravitating disks and is shown 
to yield the standard a-model in the case of non self- 
gravitating disks (Duschl et al. 2000). They have demon- 
strated that in the case of fully self-gravitating disks this 
model may reproduce very well the observed spectra of 
proto-planetary disks and yield a natural explanation for 
the accretion rate from the observed metallicity gradients 
in the disk galaxy. Furthermore, their study have pre- 
sented that the standard model of thin accretion disk 
based on a model leads to inconsistences if self-gravity 
plays an important role. This problem arises from the 
parametrization of viscosity in terms of the local sound 
speed and the vertical disk scale hight. 

Following Duschl et al. (2000) suggestion for a j3- 
prescription of viscosity, AGS06 have applied this model 
for a thin axisymmetric, polytropic, self-gravitating disk 
around a new born star. Their results is quite differ- 
ent with standard a disks in the outer part of the disks 
where the self-gravity becomes important. But, in the 
inner part of the disks their solution converged to the 
standard a disks. 

The next important possibility for angular momentum 
removal is an outflow or wind whose existence in many 
accreting systems is supported by strong observational 
evidences (e.g. Mobasher & Raine 1989; Whelan et al. 
2005; Bally 2007; Dionatos et al. 2009, 2010). Out- 
flows are generally divided into two classes, called winds 
(poorly collimated) and jets (highly collimatedfl. The 
distinction between two classes is not always clear and 
in some cases, both are seen in the same object (see Mur- 
ray 2002). It was long apparent that a disk wind/outflow 
contributes to loss of mass, angular momentum, and 
thermal energy from accretion disks (e.g. Piran 1977; 
Blandford & Payne 1982; Pudritz 1985; Konigl 1989). It 
also appears to be an almost universal feature of disk- 
accreting systems on all astrophysical scales and frame- 
works. The footprints of mass loss are observed around 

1 In this work, we will mainly focus on winds, and the study of 
rather energetic winds or jets, whose speeds are sometimes compa- 
rable to that of light, is beyond the scope of the present study. 



microquasars and YSOs as well as around massive stars 
(see chapter by Arce et al.) and even brown dwarfs (e.g. 
Ferrari 1998; Mirabel & Rodriguez 1999; Bally et al. 
2007; Whelan et al. 2005; Bourke et al. 2005), im- 
plying that the mechanism is of importance across the 
entire stellar mass spectrum (see review by Pudritz et al. 
2007). It is now widely accepted that winds or outflows 
have their origin in accretion flows (e.g. Blandford & 
Payne 1982; Fender et al. 2004). Given that the major- 
ity of mass is accreted during earlier embedded phases, 
understanding disks at early times is more critical in a 
general sense. The initial winds could be blown from col- 
lapsar disk at these early phases, when the central YSO 
still has only a fraction of a solar mass (e.g., see, Banerjee 
k Pudritz 2006, 2007). 

The accretion flows lose their mass by winds as they 
accrete onto the central object. Here, it is found that the 
mass loss rate is dependent on radius of disk and stellar 
mass, and in order to study the dynamics of the disk- 
wind system, we shall describe the shape of this depen- 
dence by a power law. In the simplest standard model 
of star formation without wind (Shu 1977; Terebey et 
al. 1984), collapse proceeds in an inside-out fashion, be- 
ginning in the center of the core, moving outward at the 
sound speed, and giving rise to a constant mass accretion 
rate of ~ 2 x 10~ 6 M Q yr _1 . Of course, many modifica- 
tions to this model have been also explored during last 
decades (Dunham & Vorobyov 2012). As a result of mass 
loss, the accretion rate is no longer constant in radius, 
and has a power low dependance on it, with the power 
low index treated as a parameter of order unity (e.g. 
Blandford & Begelman 1999; Abbassi et al. 2010). In 
the earliest phase of stellar evolution, average accretion 
rates are very high (nearly larger than 10 -6 M Q yr _1 ) 
and disk winds accompany accretion with a wind mass 
loss rate that scales as ~ xO.l the accretion rate (Gorti 
& Hollenbach 2009). Our understanding of the engine 
that powers the winds is, however, still limited. As for 
the generation mechanisms of the disk winds, there are 
several important ingredients, such as electric field gen- 
erated by the relative separation between ions and elec- 
trons, effects of magnetic fields, collisionless versus com- 
plications such as the effects of electron-positron pairs 
and coupling with the radiation field in the disk winds 
(Takahara et al. 1989). Accordingly, to date, various 
driving sources of winds are proposed, including ther- 
mal, radiative, magnetic and other ones. Traditionally 
the name of wind depends on its driving force. 

In vast majority of studies, authors readily suppose 
that the interstellar gas is isothermal and its pressure 
is proportional to the density. On account of isother- 
mal collapse, the gas had to be assumed to already be 
clumpy and at very low temperatures (a few x 10 K). 
But undeniably, the first and most powerful break from 
this isothermality comes from protostellar radiation at 
the preliminary phases of star formation (Hansen et al. 
2012). Examples include massive protostars which are 
capable of heating an entire cloud (Krumholz et al. 2007; 
Cunningham et al. 2011; Myers et al. 2011), and in 
lower level, low-mass protostars which of course may not 
have a same long-range influence as massive kinds. In the 
present paper we presume the pressure to be proportional 
to the power of density (polytropic condition) , instead of 
being proportional to the density alone (isothermal con- 



VISCOSITY AND WIND ON POLYTROPIC SELF-GRAVITATING DISKS 



3 



dition). Then, we compare the collapse of a polytropic 
gas disk with that of an isothermal one. We shall con- 
centrate on the disk structure in the presence of rotat- 
ing outflows or winds near newborn objects. However, 
since these winds are presumably launched by massive 
accretion disks (which is definitely true in early stages of 
protostellar evolution), we also pay heed to the possible 
influences of self-gravity on disk-wind system. 

Anyway, the influence of winds on accretion disks has 
been investigated by several authors (e.g. Meier 1979, 
1982; Fukue 1989; Takahara et. al 1989; Knigge 1999, 
hereafter CK99; Xie & Yuan, 2008; MS09; Kawabata 
& Mineshige 2009; Dotan & Shaviv 2011; Abbassi et 
al. 2010). In these studies, the accretion disk is usually 
height-integrated and the configuration of the accretion 
flow is assumed rather than being calculated. Jiao & Wu 
(2011) solved the full hydrodynamic equations to get the 
configuration of the accretion flow. CK99 derived the 
radial distribution of dissipation rate and effective tem- 
perature across a Keplerian, steady-state, mass losing ac- 
cretion disk, using a simple parametric approach. MS09 
studied the influence of the winds on the time evolution 
of isothermal, self-gravitating accretion disks by adopt- 
ing a power law mass loss rate because of the existence 
of wind. The work by MS09 considered both mass and 
angular momentum loss due to the wind/outflow, and 
provides a basis for the present work. 

In accord with earlier discussion, viscosity, self-gravity 
and also outflows are among the most important phys- 
ical agents in accretion-dominated stages of star forma- 
tion. Although a considerable amount of work has gone 
into any pair of these, but no study has combined all 
three with non-isothermal circumstances using similar- 
ity approach. To our best knowledge, the simultane- 
ous solutions for a viscous (with two prescriptions) , self- 
gravitating, polytropic, mass-losing disk have not yet 
been reported. 

Now, we are interested in considering the possibil- 
ity that winds, could affect the global properties of 
polytropic self-gravitating accretion disks. Indeed, the 
present work is an attempt to provide a more thorough 
survey of solutions for two viscosity models with poly- 
tropic condition and wind, adopting different values of 
input parameters. As CK99 and MS09 the parametric 
model adopted here to describe the mass loss is simple, 
yet sufficiently general to be applicable to many types 
of dynamical disk-plus-wind models. The paper orga- 
nized as follows. In the next section we describe a basic 
physical approach for setting up analytical part of our 
model. We did not consider the driving mechanisms of 
the wind. In §3 we firstly attempt to reduce our similar- 
ity equations in a slow accretion limit. Then, assigning 
both proposed models for viscosity, we find two sets of or- 
dinary differential equations and explore their solutions 
in §4. A summery on the properties of our solutions and 
their implications are given in §5. 

2. GENERAL FORMULATION 

2.1. Basic equations 

Let us suppose a geometrically thin accretion disk sur- 
rounding a central object which has not yet been com- 
pletely formed. Thus, the radial component of the gravi- 
tational force is mainly provided by the self-gravity of the 



disk itself. The disk we are concerned with in this paper, 
has a symmetry over the rotation axis and is turbulent 
having an effective turbulent viscosity. Besides above- 
mentioned effects, we consider a wind emanating from 
disk surface, and set up fundamental governing equa- 
tions in the cylindrical polar coordinates (r, tp, z) cen- 
tered on accreting object with the equatorial plane of 
disk at z = 0. The equations of continuity, radial mo- 
mentum and conservation of angular momentum, respec- 
tively read (e.g. CK99; MS09) 



do | 1 d | 1 

dt r dr r 2irr dr 



0, 



~dt 



dv v 
dr 



r 



ldp 
p dr 



(1) 

(2) 



d(r 



dr 



ld_/ 3 dQ.\ (lr) 2 n dM* 
r dr V dr J 2irr dr 



d(rv v ) 

(3) 

where v T , v v , fl, p, o, p and v are the radial, rotational 
and angular velocity, volume density, surface density, gas 
pressure and kinematic viscosity of the disk, in the or- 
der given. Other variables are set out in detail shortly 
afterwards. 

Insofar as we require all quantities of the flow vari- 
ables depend only on radius and time, we have integrated 
them in the vertical direction, following the thin disk ap- 
proximation. So, rather than dealing with quantities per 
unit volume (such as the density p) , we deal instead with 
quantities per unit surface (such as the surface density 

In Eq. ([2]), we have radial component of the gravita- 
tional acceleration due to the self-gravity (cf. Saigo & 
Hanawa 1998; TT1999), 



<9$ 

dr 



OO />OQ 



= -2ttG / / J 1 (kr)J (kr')o(r')kdkdr' , (4) 
Jo Jo 

derived from the axisymmetric Poisson equation adapted 
for a thin disk in cylindrical space (Nomura & Mineshige 
2000), 



IS / 9$ 



r dr\ dr J 



<9 2 $ 

dz 2 



4irGo(r)d(z). 



(5) 



Here, functions J\ and Jo are the Bessel functions of 
the first kind and 5(z) is the Dirac delta function. As 
usual, G is the universal gravitational constant and $ is 
the gravitational potential contributed by the material 
of mass 

(6) 



M r (r) = 2tt / cr(r') r'dr', 
Jo 



inside of the disk within cylindrical radius r. 

Since the gaseous disk is not rotating as a solid body, 
in the azimuthal equation of motion [Eq. (j^J], we have 
the shear viscosity whose presence allows the transport 
of angular momentum from the faster inner fluid ele- 
ments to the slower outer ones. Regardless of mecha- 
nisms which might contribute to the initiation of an out- 
flow from the surface of an accretion disk, mass loss rate 
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by wind/outflow is represented by M w in Eqs. ([TJ and 
([3]), and it is newly included in our analysis, whereas the 
Euler Eq. ((2]) does not change anymore despite inclu- 
sion of the wind. It is easy to perceive that merely the 
vertical component of wind velocity would make the gas 
escape from a surface of the disk. In fact, we have 



M w (r) = 2 



r (r') 2-Kr'dr', 



(7) 



where a w — pv+ is mass loss rate per unit area from 
each disk face, picking p = p| z= o as a midplane density 
of the disk and v£(= —v~) > as mean vertical veloc- 
ity at the disk surface, i.e. at the base of a wind. A 
factor of 2, appeared behind the integral, is for taking- 
account of ejection from both disk faces. Determining 
the wind velocity at its base is complicated because it 
depends on the vertical structure of the disk which in 
turn depends on the detailed variation of the unknown 
viscosity with height. Furthermore, the opening angle for 
the wind flow needs to be determined using the geome- 
try of the accreting flow and the pressure gradient in the 
wind in the radial direction (Misra & Taam 2001). Be- 
cause is unbeknown to us, we will compensate this by 
constructing a library of wind solutions for a wide range 
of wind model parameters with v+ entangled therein. 

The rightmost term of Eq. ((3]) is the outflow sink added 
term represents angular momentum transferred by the 
wind. Here, it is assumed that matter ejected at radius r 
on the disk, carries away the angular momentum (Zr) 2 f2, 
where Q = v v /r is the angular frequency associated with 
lever-arm, Ir, at radius r. Thus, I = corresponds to a 
non-rotating disk wind and I = 1 to outflowing material 
that carries away the specific angular momentum (r 2 fi) 
it had right at the point of ejection. This latter would be 
the most fitting value for radiation-driven outflows (e.g. 
Murray & Chiang 1996; Proga et al. 1998; Feldmeier 
& Shlosman 1999; Feldmeier et al. 1999). Moreover, 
< Z < 1 hints at the family of disk winds that carry 
away less angular momentum than possessed by the wind 
material before it left the disk surface. Centrifugally- 
driven MHD disk winds (magnetocentrifugal winds for 
short) are corresponding to / > 1 and can remove a lot 
of angular momentum from the disk (e.g. Blandford & 
Payne 1982; Cannizzo & Pudritz 1988; Emmering et al. 
1992; Pelletier & Pudritz 1992; Pudritz et al. 2007). In 
this case, we have I = Tx/t, where ta is Alfven radius 
(CK99). The angular momentum that is observed to be 
carried by rotating flows (e.g. DG Tau) is a consistent 
fraction of the excess disk angular momentum, from 60- 
100% (e.g. Bacciotti 2004), due to the high extraction 
efficiency mentioned above. As we will discuss later in 
§5, it is observationally well-known result that in many 
systems, M w /M acc ~ 0.1, with M acc being mass accre- 
tion rate which is introduced in £ )2.2I This is faithful to 
the fact that lever-arm coefficients under the latter pic- 
ture are often found in numerical and theoretical works 
to be / = ta jr ~ 3 - the observations of DG Tau being 
a perfect example (see a review by Pudritz et al. 2007). 

To concede the truth of temperature distribution dur- 
ing the collapse, we employ a polytropic relation between 
pressure and density of accreting gas 



p = Kp 1+ » = Kp~< 



(8) 



where A" is a constant set by the entropy of the gas and n 
is known as the polytropic index for the process of inter- 
est. It is customary to replace 7 by 1 + i and call it an 
effective adiabatic index. We should imply here that in 
actual accretion flows that can be represented by poly- 
tropes, neither K, nor n, may not be constant in space 
or time, and j(= dlap/dhip) does not necessarily equal 
1 + i (e.g. Goldreich & Weber 1980; Yahil 1983). But we 
shall confine attention, in this paper, to an abstract case 
where each has certain fixed amount. The underlying 
polytropic exponent 7 is somehow an adjustment screw 
by which we can tune how much the temperature varies 
throughout the flow. The temperature increases more 
rapidly towards the center of polytropic gaseous disks 
with larger values of 7 (see MNU97). As 7 approaches 
unity, the temperature inclines to be uniform all around 
the disk, and consequently all values in the disk tend to 
those of isothermal case. The values in the disk such as 
geometrical thickness depends critically on 7 (see Saigo 
et al. 2000). According to Omukai & Nishi (1998), the 
polytropic relation with 7 1.1 is a good approximation 
to a collapsing metal-free gas cloud (see also Matsuda et 
al. 1969; Carlberg 1981). If c s is the local mean sound 
speed in barotropic disk, it satisfies 



2 d P is 7-1 
c s = — = K^p< . 

dp 



(9) 



Using last two equations, the pressure force per unit mass 
on the right-hand side of Eq. ^ becomes 



1 dp c 2 dp 
p dr p dr 



a dr 



(10) 



From the hydrostatic equilibrium equation in the vertical 
direction, 

c 2 da d<$> 



a dz 



dz 



0. 



11 



and the Poisson Eq. ([5]), we can get the vertical extent 
of the disk at any radius, as 



H 



(4ttGp)2 



27rGo-' 



(12) 



viz. half-thickness of the disk, where we assume that the 
differential of the gravitational potential in the radial 
direction is negligible compared with that in the vertical 
direction (e.g. Nomura & Mineshige 2000). Here, as 
stated, the azimuthal integration leads to 



o-(r) = 



+H 



p{r,z) dz « 2p(r,0)H(r). 



(13) 



H 



At last in this subsection, we apply the latest changes 
in our fundamental Eqs. (JTJ)- (J3J) and rewrite them as 



da 1 d 



dt 



2i 

dt 



dv T 
~dt 

dj 



i'r 



r dr 

dv r 
dr 



(rav T ) + 2a v 



0, 



c[da_ 
a dr 



3r, 



(14) 
(15) 



_ J_d_ 

r dr ra dr 



- -a 2 ^w), (i6) 

or r z 
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where the angular momentum j is replaced by rv v , and 
we used Eq. ([7]) to obtain the radial derivative of mass 
loss rate, i.e. our auxiliary equation 



dM* 
dr 



(17) 



The solutions of recent four equations give us an in- 
clusive overview of the evolution process in infant stellar 
objects, which strongly depends on the chosen viscosity 
model, and partly on the strength of existing wind. So, 
the study of dynamical behavior of the accretion disks is 
postponed to more information about the viscosity and 
also mechanisms through which the mass and angular 
momentum loses to the wind. 

2.2. Self-similar scaling 
The basic Eqs. (| 14[) - {(T6|) are troublesome to solve in 



real space, albeit approximately. For a sense of how they 
might be solved, it is useful to alter them to their non- 
dimensional forms, that firstly involves adopting an ap- 
propriate similarity variable. In this way, with the aid of 
K and G as the constitutive dimensional parameters, we 
pose a dimensionless similarity variable x = f sc r, with 



f sc = K- L 2G i —(±tr- 2 



(18) 



being the inverse length dimensional scale factor (e.g. 
MNU97, Yahil 1983, AGS06). Using such a combina- 
tion of radius r and time t, all hydrodynamic variables 
must therefore be functions of x only and solutions will 
have the spatial structure at all times, because of self- 
similarity. Once the core formation epoch is an origin 
of time (t — 0), we take the plus sign in the parenthe- 
sis so as to gain positive time (t > 0) solutions. Hence, 
bearing in mind the chain rule for the transformation 
(r,t) —> {x,t'), derivatives will turn into 



or dx 

d d , .x 8 
77 + (7- 



dt dt' 



f dx 



(19) 



(20) 



From now on, we are allowed to write t instead of t' 
and d/ dx instead of d/dx, since we require that all time- 
dependent terms should disappear in the self-similar 
forms of equations. Now, it is convenient to draw up 
a list of useful real quantities correlated with their self- 
similar kinds: 



p{r,t) = {A^y-G^r^-ix) 
p(r,t) = {^ 1 )- 1 KG-^r 2 ^ 2 {x) 
c s (r,t) = (4 7 r)^7^idG i ^t 1 - 7 £ J T i (2) 
v r (r,t) = K^G^t^VAx) 
v v {r 1 t) = K^G^ L t 1 -^V v {x) 
v+(r,t) = K^G^t^V+ix) 
j(r,t)=](G H ^ 27 J(x) 
<r(r,t) = (2n)- 1 KiG- 1 ^ L t^j:(x) 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



H{r,t) = ^tOV^^sgV^^V^) (29) 



v{r,t) =KG 1 -"<t 3 - 2 " / v'(x) 



(30) 



M T (r,t) = KiG^t^M^x) (31) 



M r (r, t) = KiG^t^-^M^x) (32) 



M W M) = KiG^t^-^M^x) (33) 



dw(r,t) = (4ir)- 1 K^G- 1 ¥t-<- 1 Z(x)r(x) (34) 



where 



J(x) = xV v (x), 
M K (x) = [ E(x')x'dx', 



M w {x) = / T,(x')T(x')x'dx' 
Jo 



(35) 
(36) 

(37) 



T(x) = (47r) 1: r L 7 -^E Hi V li ( a ;)A(x) , A = EV+. (38) 

Also, the velocity of the constant x surface seen from the 
rest frame (r, t) is 



Urest(V,t) 



dr 
dt 



V Icst {x) = {2- 1 )x. (39) 
For later convenience, we define the comoving velocity as 

V = V T - Kest =V t + (-y- 2)X. (40) 

It is useful to rewrite mass conservation Eq. (TT]) in terms 
of enclosed mass M r [Eq. ((BJ], i.e. 

8M T dM r ■ n dM T „ 

+ v r ^- + M w = , = 2-nra. (41) 



dt dr 



dr 



Then, introducing M acc as the mass accretion rate, we 
have M acc = M T + M w , in such a way that 



M acc = -2irrav r . 



(42) 



A minus sign in this expression obviously shows when 
the radial velocity of the gas flow is directed inwards, an 
accretion may take place. Self-similarity leads to 

M acG (r,t) = KiG^t^-^M^dx), (43) 

ExV v 



Me 



37-4' 



(44) 
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2.3. Basic equations in self-similar space 

After some replacements with the contribution of pre- 
vious subsection, it is now possible to recast the differ- 
ential Eqs. (jT^T| - ([TT)) . from partial (PDEs) into ordinary 
ones (ODEs). The nondimensional similarity equations 
are then derived as 



— (EVx) + (r - 3 7 + 4)Ea: = 0, (45) 
dx 



^ + ^-4-<?x-(2 7 -3)V+(7-2)( 7 -l)z = 0, 
ax 1j ax x a 



(46) 



V 



dJ Id 

dx Ex dx 



z/a: 3 E 



d i J 



dx \x 

dM v 
dx 



(^)]+(' 2 r-2 7 + 3)J = o, 

(47) 



Ear, 



(48) 



where, for conciseness of notation, we introduced b = 

1 y 2(7-1) 

(47r 7 ) i E i , and the similarity gravitational field in 
radial direction is 



(49) 



G x (x) = - / / J 1 (kx)J (kx')T,(x')kx'dkdx'. 
Jo Jo 

Clearly, the effect of wind or outflows appears by the 
term T. When we set this parameter equal to zero, the 
fourth equation can be left out and the other ones reduce 
to those appeared in MNU97 for a-case or AGS06 for (3- 
case (see 93 ,2[) . Although full numerical solutions to these 
equations would now be possible, it is more instructive to 
proceed by analyzing the model in some restrictive cases 
such as one on the slow accretion limit. 

3. REDUCTION OF BASIC EQUATIONS 
3.1. Helpful approximations 

In this subsection, we lay eyes on two useful approxi- 
mations, and see how they reduce and simplify our equa- 
tions, without much loss of generality of the problem. 

Firstly, it seems favorable to ignore some terms like 
pressure gradient force and acceleration term in the cold 
and slow accretion limit which implies ^ > 1, E > 1 
and \V r \ -C 1. This estimate has been widely used by 
many authors to simplify formulas (e.g. MU97, MNU97, 
TT99, AGS06, MS09). Furthermore, because in our 
model the wind velocity is expected to be smaller than 
radial inflow velocity, one can handily assume M w < 1 
as another implication of this limit, whose consistency 
with the results can be readily verified by reader. We 
will also check that M w do not exceed the accretion rate, 
M acc . From a different but equivalent point of view, the 
slow accretion approximation is applicable for rotation- 
ally supported disks when the viscous timescale is much 
longer than the dynamical timescale (TT99; AGS06). 

Next, to avoid the integro-differential equation arises 
from substitution of Eq. (|4T?)) into Eq. (|4"6")l . we use a 
monopole approximation to compute the radial gravity 
field consistent with the mass distribution. Adopting this 



approximation under the restrictions of the slow accre- 
tion limit, leaves (cf. TT99) 



-EV 



(3 7 - 4)x 



(50) 



Except near the outer edge, this neglect of higher mul- 
tipole orders is not expected to introduce any signifi- 
cant error as long as the surface density profile is steeper 
than l/r (e.g. MU97; MNU97; Li & Shu 1997; Saigo 
& Hanawa 1998; TT99; Krasnopolsky & Konigl 2002; 
AGS06; MS09). 

In the slow accretion limit, the third and fourth terms 
in Eq. (|4"6")l dominate and the others could be canceled 
(e.g. MNU97; AGS06). Then, by making use of last rela- 
tion, we can demonstrate a radial force balance supplied 
by two terms in Eq. (|46|) . 



J 



x 3 (3 7 — A)x 



0, 



leading to 



J 



(—)' 

\3 7 -4/ 



(51) 



(52) 



One can take a logarithmic derivative of J with respect 
to x, and get 



din J 
dhix 



ld\n\V\ 



2 d\nx 2 d\nx 



1 dlnE 
1 + t: 



(53) 



After some algebraic manipulations, Eq. (|4"5"j) takes the 
form 



rflnE _ d\n\V\ 



(37- T- 4)-, (54) 



d In x d In x 

then can be substituted in Eq. (f53")l . to obtain 

din J 1 . /3 7 -T-4 



din a; 



2V 



(55) 



To proceed, from this point on, the kinematic coeffi- 
cient of viscosity v' needs to be assigned as a function of 
our similarity variable x. 

3.2. Viscosity prescription 

For the viscosity v' we shall apply two different rep- 
resentations. The first and obvious one is the a- 
prescription introduced by Shakura & Sunyaev (1973) 
and the resulting polytropic a-disk plus wind will be for- 
mulated in the next subsection. As mentioned in JT] the 
a- viscosity is not a unique choice and we also employ the 
so-called /3-viscosity introduced by Duschl, Strittmattcr 
& Biermann (2000) thereinafter. 

3.2.1. ce-model solution 

In the case of a-prescription, as suggested by MNU97 
we adopt v' a = a'x , with a' and I being free parameters. 
So, the viscosity coefficient is a function of x only. Sub- 
stituting this prescription into Eq. (|47j) we can inquire 
into a dynamical evolution of the disk. 

Eqs. (|S1|) and ([5SJ help us to simplify Eq. (gTJ) and 
after some mathematical manipulations we can finally 
obtain a desired first order ODE for a-viscous disks as 
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dV _ 1 [K- + (2Z 2 -l)aT][K- + (7-2M 
dx a ~ u' a 3F r - (37 - 2)x + 2xT 
Ax 2 + BV T x + 3(21 - 1)V 2 
+ 2x[W T - (3 7 - 2)x + 2xT] ' 

where for a sake of brevity we have used 



(56) 



A = A{2i + 3) - A{1 + 3)7 - 3r 2 
+ [2£(7-2) + 6 7 ]r-2(2- 7 ). 



dT 

dx ' 



B = 6(£ + 2)7 - 8(2£ + 1) + 2{l - 4)r + 2.x 



dT 

dx 



Clearly, the effect of wind/outflow appears by the term 
r. If we set this parameter equal to zero, we can find the 
Eq. (34) of MNU97 which matches the solution without 
wind. Considering £ — 1 and 7 = 1 for an isothermal 
case, the equation reduces to Eq. (19) of MS09. If wc 
continue by setting r = 0, we could easily obtain the Eq. 
(18) of MU97. 

3.2.2. /3-model solution 

As we discussed in introduction we are also willing to 
use so-called /3-prescription as an alternative model for a 
self-gravitating disk which is introduced by Duschl et al. 
(2000) and is used by AGS06 in the form 



y /3 



p'xV v = P'J, 



(57) 



with j3' being a free parameter. It is worth noting that 
Duschl et al. (2000) /3-viscosity recovers Shakura & Sun- 
yaev (1973) cn-viscosity for non-selfgravitating disks, if 
one requires the turbulence not to be super-sonic (which 
makes sense and is in agreement with Shakura and Sun- 
yaev's parametrization, but still, it is an additional con- 
dition which one should mention). Likewise, with the 
aid of Eq. ([55)1 and last equation, from Eq. (|47p we can 
derive the single first order ODE for /3-viscous disks as 



dVr 
dx 



P 



1 [V r + {2l 2 - l)xT][V r + (7 - 2)x] 
t/p 3V r - (37 - 2)x + 2xT 
CV r + Dx 



W r -(3-/-2)x + 2xT~ 
where again for short we have written 



(58) 



C = 3(2 - 7) + 5(37 - T - 4) + x 



dT 

dx ' 



D=(2-j)(7T-x 



dT 

dx 



I87 + 22) - 2(3 7 - 4) 2 



Same as previous, one may ignore the influence of wind 
in /3-model, by adopting T — 0, and this time recover the 
Eq. (32) of AGS06. 

The ordinary differential Eqs. (1551) and ([55} along with 
Eq. (j48|) and Eq. (|54[) are the main equations of our 
analysis which we can solve them numerically using the 
forth-fifth order Runge-Kutta-Fehlberg scheme. By ex- 
ploit of asymptotic solutions as boundary conditions for 
our equations, we will investigate the effects of physical 
parameters on structure of the disk-wind system. 



4. NUMERICAL ANALYSIS 

In our model, there is a set of the input parameters like 
s, I, Ao, 7, a' and /3'. Thus, a clear physical picture of our 
model is obtained only by an extensive parameter study. 
Here, we first obtain appropriate boundary conditions, 
and the possible effects of our different input parameters 
are explored in the subsequent subsections. 

4.1. Boundary conditions 

Our derived ODEs require boundary conditions so as 
to be solved numerically. Thus, at the first step, we 
derive asymptotic solutions for V r . Our first limit is near 
the origin of the disk where x 0, and the second one 
is in the outer part of the disk that, i.e. x — ¥ 00, which 
is toward the parent cloud. As a result for a-model we 
have 



lim Vr 



lim V r a ' 

x— > + oo ' 



(4£ + 6) - {21 + 6)7 
2,-21 

9(4- 3 7 )x 3 



(2l--S)(M 2 -2£-\h)v' a 
AvL 



2x(j - 2) 2 
Also for /3-model it can be written as 



(2l J - l)xT 



lim Vr g ' 



lim Vr 3 ' 

x— >+oo ' 



57 - 6 r (7 - l) 2 a; 3 



67-7 
x{ 1 -2f 



(67 - iy B 

(2l 2 - X)xY 



(59) 
(60) 

(61) 
(62) 



As 7 approaches unity, the results approach the isother- 
mal case, again as expected. Although to start the in- 
tegration one can assume V r — at x — 0, we should 
avoid the singularity by avoiding the origin. The above 
asymptotic behaviors helps us to determine an appro- 
priate boundary conditions at a small neighborhood of 
singular point x = 0. To have surface density profile, we 
can integrate surface density equation [Eq. (|54|) ] from 
outer boundary, i.e. x = 1, towards the center of the 
disk for a given 'Six — 1) = E out . Besides this, we de- 
mand a series of modes in which the accretion rate at the 
outer edge (i.e., the inflow rate from the parent cloud) is 
kept constant, which appears as a natural requirement 
(MS09). Thus, our main boundary conditions are V r = 
at x = 0, and, M acc = M infaU at x = 1. So, Af in f an is an- 
other input parameter to be assigned to our model. For 
the profile of mass loss rate, we prescribe a basic power 
law form for A as A = Aqx s , with Aq and s being the free 
parameters. Note that in the case of isothermal collapse, 
one derives A = T, and the wind parameters reduce to 
those of MS09. 

Now applying boundary conditions for our two systems 
of ODEs, we can obtain some profiles for hydrodynamic 
variables. Figs. [T][5] show radial distributions for some 
important physical variables with parameterized values 
as a function of similarity variable x. For instance, one 
would immediately deduce that the radial velocity in- 
creases meaningfully at the outer part of the disk, be- 
cause of the wind. In order to make an easier comparison 
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between surface densities, the ratio (a — cro)/co versus x 
is shown in all figures and its negative value means that 
the surface density generally reduces in the presence of 
the disk wind. In the figures, we see the mass accretion 
rate plainly decreases in comparison to the no-wind so- 
lution. It is also an easy task to compare a with /3 disks. 
In the outer part of the disk - where the self-gravity is in- 
fluential - the behavior of the solutions predicted by the 
f3 viscosity model shows much less radial velocity com- 
pared to that by a viscosity model, either for wind or 
no-wind case. Moreover, it is apparent from the figures 
that winds could lessen the rotational velocity of the disk. 
Wind solutions imply that the amount of reduction to the 
rotational velocity is more significant for /3-disks than it 
is for a-disks, under our restrictive boundary conditions. 

4.2. Role of mass loss index s 

Figs. [T]andE]show the effect of adopting various values 
for the mass loss power law index s, on the profiles of the 
physical variables. Each curve is labeled by correspond- 
ing index s. Also, we adopt a' = 0.1, £ = 1, /3' = 10~ 3 , 
M infaU = 2.0 x 10~ 6 M Q /yr, A = 0.1 and I = 1 (i.e., 
rotating wind). Our adopted value for M; n f a n is also 
compatible with the mean values inferred in embedded 
protostars (e.g. Konigl & Pudritz 2000). The surface 
density and the rotational velocity for the no-wind solu- 
tion are represented by ao and vq, v , respectively. In order 
to make an easier comparison, the ratio (a — qo)/ 'o-q as a 
function of variable x is shown in Figs. Q] and [5] (middle, 
left). Since the similarity radius is smaller than unity, for 
smaller values of parameter s (consider A = Aox s ), the 
wind becomes stronger and more mass is extracted from 
the disk. Therefore, surface density reduction is more 
significant for smaller values of s. 

The rotational velocity profiles are presented as a func- 
tion of the similarity variable in Figs. [1] and [2] (top, mid- 
dle). Generally, when the exponent s decreases, the flow 
will rotate slower than that without winds. So, the vis- 
cous dissipation per unit mass in the flow is expected to 
be smaller in the presence of a wind. Also, the radial 
velocity profiles in Figs. [T] and [5] (top, left) represent sig- 
nificant deviations from no-wind solution. In the outer 
parts of a disk the radial velocity is approximately uni- 
form in the no-wind case. But, as the wind plays its 
crucial role at the outer parts of the disk, we have much 
larger radial velocity in comparison to the no-wind so- 
lution. Wind velocity w+ at the surface of the disk is 
shown in Figs. [1] and [5] (top, right). /3-disks considerably 
have less v£ than a-disks. As wind gets stronger, its ve- 
locity at the surface of the disk increases which is quite 
expectable. 

The accretion rate profiles, M acc , are shown in Figs. [1] 
and [2] (middle, middle). Accretion rate for the no- wind 
solution is represented by the dashed curves. Generally, 
the accretion rate decreases at all parts of the disk in the 
presence of the wind. Nevertheless, the accretion rate is 
not much sensitive to the variations of the exponent s. 
Ratio of the mass loss rate by wind to the accretion rate, 
i.e. the mass loss efficiency M w /M acc , is plotted in Figs. 
Q]and[5] (middle, right). The mass loss due to the wind 
is negligible in the innermost region of the disk, except 
for strong winds, e.g. for s — 0.1. Our input parameters 
were chosen so that mass loss efficiency is less than one at 



all radii of the disk. Here, the larger mass loss efficiencies 
appear at large radii, i.e. outer part of the disk. For a 
stronger winds which correspond to smaller values of s, 
a larger fraction of the mass carries away by the wind. 

The profiles of disk aspect ratio, H/r, for various val- 
ues of exponent s are presented in Figs. [1] and [2] (bot- 
tom, left). They demonstrate that in general a-disks are 
thicker than /3-disks and both of them are thicker in the 
presence of wind with all adopted values for s. In Figs. 
Q] and [2] (bottom, middle), we can see the distribution of 
angular momentum per unit mass, i.e. j/M T = j, along 
the similarity radius of the disk, x, for different expo- 
nents s. As winds emanate from the disk, this fraction 
would be larger, particulary for smaller values of s. One 
would see the same behavior in a /3-disk that typically 
has the lower values. Although there is an angular mo- 
mentum loss due to the wind, the accompanying mass 
loss is high enough to keep j increasing as wind becomes 
stronger. 

4.3. Role of dimensionless lever-arm I 

Possible effects of an angular momentum extraction 
due to the wind are explored in Figs. [3]and|4]by adopting 
various values of the input parameter I. We here assume 
a' = 0.1, t = 1.0, p = 10" 3 , M infaU = 2.0 x 10- 6 M Q /yr 
and I = 0.0, 1.0, 1.5 with A = 0.1 and s = 0.7. Ob- 
viously, when we have 1 = 0, angular momentum is not 
extracted by the wind. This case corresponds to a non- 
rotating wind and the disk losses only mass because of 
the wind. However, as mentioned in MS09, it can be 
shown that for I 2 < 1/2 the mass of the disk increases 
in the presence of the winds that obviously has not any 
physical meaning. This is partly due to the limitations 
of similarity method that there is not a self-consistent 
solution for any given set of the input parameters. More 
importantly, our model is valid just in the slow accretion 
limit which implies V v ^> 1, and so it is very unlikely to 
accept that winds are lunched without extracting a cer- 
tain amount of angular momentum of the disk (MS09). 
Although we have represented solutions with I = 0, 0.5 
in Figs. [3] and [4] for a sake of comparison, as in MS09, 
we think these solutions are not physically acceptable. 
Profiles of surface densities for each viscosity model are 
shown in Figs. [3] and 0] (middle, left). We can see again 
the reduction to the surface density because of the wind. 
Rotational velocity of the disk decreases because of the 
angular momentum removal from the disk, as I becomes 
larger (top, middle in figures). Although, the radial ve- 
locity in the innermost part the disk does not change 
because of the wind, but in comparison to the no-wind 
solution, existence of a rotating wind enhances the ra- 
dial velocity at the outer part of the disk (top, left). The 
typical behavior of wind velocity v% is also sensitive to 
the amount of the extracted angular momentum (top, 
right). The accretion rate profile (middle, middle) repre- 
sents that it decreases due to the existence of a rotating 
wind. However, for a large I, where more angular momen- 
tum is carried away by the wind, as long as the surface 
density and the rotational velocity are reduced at all re- 
gions of the disk, the radial velocity of the accretion flow 
at the outer part of the disk is increased significantly 
(cf. MS09). We also plotted the specific disk angular 
momentum j as a function of x so that one can see its 
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behavior for a and (3 disks assuming 7 = 1.1 for different 
values of I. Here, it is informative to compare solutions 
for I = 0, 0.5 with those for I = 1.0, 1.5. It can be inferred 
from the disk aspect ratio profile (bottom, left) that the 
more angular momentum is removed from the disk, the 
thicker it gets. 

4.4. Role of factor A 

One of the prominent input parameters in our model 
is Ao that its possible effects arc explored in Figs. [5] 
and El We assume that a' = 0.1, I = 1.0, /3' = 10~ 3 , 
M infaU = 2.0 x 10- 6 A/ Q /yr and A = 0.1, 0.05, 0.01 with 
I = 1.0 and s — 0.7. The surface density and the rota- 
tional and radial velocities are substantially decreasing 
with Ao. The wind velocity at the surface of the disk is 
significantly affected by the parameter Ao (top, right). 
As a result, the mass accretion rate and the wind mass 
loss rate are respectively decreased and increased with 
the parameter Ao. Further, j increases in both viscos- 
ity models, as wind becomes stronger by adopting larger 
values of A . 

4.5. Effect of disk self- gravity 

The effect of the disk self-gravity in this paper is lim- 
ited to provide the radial gravitational field to keep the 
disk in centrifugal equilibrium. On the other hand, it is 
predicted that in the outer part of accretion disks around 
QSOs, self-gravity has a dominant role. This effect is in- 
vestigated by Toomre (1964). As a simplest indicator 
for gravitational stability of the solutions we can use the 
Toomre criterion, 



where 

K = 0(4 + 2^)4 (64) 
dlogr 

is the epicyclic frequency at which a fluid element oscil- 
lates when perturbed from circular motion. In a nearly 
Keplerian disk, k ~ f2. For axisymmetric disturbances, 
disks are stable against the gravitational fragmentation 
when Q > 1. The local gravitational instability occurs 
when Q < 1. Now, we rewrite our gravitational instabil- 
ity parameter in the self-similar form as 

Q = 2V2(A7r)^ 1 ^^x- 2 j(^-Y, (65) 

Vain xJ 

which by setting 7 = 1, we recover Eq. (18) of TT99, viz. 

Q = (2V2Jy/d In J/d In x) /Zx 2 . In all fi gures (bottom, 
right), we have shown the distribution of the Toomre 
Q- value for some parameters. 

To make an easier comparison, Toomre parameter for 
a case without wind/outflow is also represented. The so- 
lutions indicate that Toomre parameter increases with 
winds or outflows, except for the cases with I = 0, 0.5 
which gives unphysical solutions, as we discussed previ- 
ously. Generally in a-model, except for the inner part of 
the disk, Toomre parameter is still larger than one, espe- 
cially when the winds are present. However, in the case 
of /3-prescription for viscosity, we see that Q is below the 
instability threshold (i.e. Qthr ~ 1) m most regions of 
the disk, even in the presence of a typical wind (but some 
cases, e.g. for I = 1.5). Manifestly, one will not have any 



trouble adjusting the input parameters for /3-disk to be 
locally unstable in various distances from the central ac- 
cretor and so the /3-disk model might be a good nominee 
for the origin of planetary systems (e.g. AGS06). 

4.6. Role of polytropic exponent 7 

Another important input parameter of our model is 
7 whose possible effects are examined in Figs. [7] and 
B We assume that a' = 0.1, I = 1.0, /?' = 10~ 3 , 
Minfeu = 2.0 x 10- 6 M Q /yr, A = 0.1 and 7 = 1.0, 1.1, 
1.2 with I = 1.0 and s — 0.7. As it is shown there, the ra- 
dial velocity, surface density, vertical wind velocity, mass 
accretion rate and mass loss efficiency are significantly 
decreasing with 7 in both viscosity models. There is 
also a reduction to the rotational velocity in the pres- 
ence of wind, and as 7 increases in a-model or decreases 
in /3-model, this reduction seems to be greater, under 
our restrictive boundary conditions. The disk aspect ra- 
tio H/r, angular momentum per unit mass j, and also 
the Toomre instability parameter Q (introduced in ^4.5[) 
are highly affected by the given parameter 7, and con- 
siderably increase with it. Thus, aside from the adopted 
viscosity prescription, the disks with larger 7s are thicker 
and gravitationally more stable. 

5. SUMMARY AND OUTLOOK 

In the present study we have examined the influences 
of hydrodynamical winds from a geometrically thin disk 
rotating around the central object, taking account of the 
self-gravitational field of the disk gas with a and (3- model 
for its viscosity. We used the self-similar method to ob- 
tain the dimensionless form of the fluid equations, and 
then reduced them in the slow accretion limit. In order 
to describe the evolution of our disk, we derived two sets 
of ODEs for two available models of viscosity. We solved 
them numerically, by exploit of natural requirements and 
asymptotic solutions near the origin and near the outer 
edge, as the boundary conditions. Of course, it is im- 
portant to keep in mind, we had the limitation to select 
parameter 7 for essence of differential equations and the 
fact that we seek just physical solutions (see MNU97 for 
constraints). 

Some fraction of the accreted material and their an- 
gular momentum can be carried away by the wind. At 
the inner part of the disk, wind does not alter consider- 
ably the dynamical behaviors of the disk. However, at 
the outer part, where the wind appears more efficient, 
all physical variables would be much modified by wind. 
We presented the ratio of the total mass loss rate by 
wind to the mass accretion rate at each radius of the disk 
in all figures. There are some observations evidences of 
different systems which show that this ratio is around 
0.1 (e.g., Konigl & Pudritz 2000), which is in agreement 
with our results. Additionally, all the figures show con- 
sistency with our assumptions due to the slow accretion 
limit (§3.1). 

In spite of simplicity of our model in treating the wind 
and the disk itself, we think the presented semi-analytical 
results give us a better understanding of such a compli- 
cated system. Basically we had three main input param- 
eters to control the physical properties of the wind in 
a phenomenological way, i.e. s, I and Ao, and another 
input parameter related to the thermodynamics of the 
disk, i.e. 7. We did an extensive parameter study for a 
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wide range of the input parameters and the main results 
are summarized as follows: 

1. Radial dependence of the mass loss by wind was 
prescribed by a power law with exponent s. As this pro- 
file of mass loss becomes steeper with the radius, the 
accretion velocity is enhanced in particular at the outer 
layers of the disk. Additionally, the radial velocity in- 
creases with A . Naturally, in both viscosity models, if 
we keep all the input parameters fixed and decrease s 
(because < x < 1) or increase A , then more mass is 
extracted from the disk by the wind. It means that more 
angular momentum is extracted by the wind in addition 
to the turbulent viscosity which implies a larger radial 
velocity as the solutions clearly show this behavior. De- 
pending on the wind mechanism, value of I is adopted in 
our model. Larger I implies more efficient angular mo- 
mentum extraction by the wind which leads to a more 
stable disk with larger radial velocity. 

2. As the wind becomes stronger, the disk losses more 
mass and so, one should normally expect a reduction to 
the surface density of the disk. Consistent with this phys- 
ical expectation, we showed that in the presence of wind, 
surface density decreases by decreasing s or increasing / 
or Ao in both viscosity prescriptions. 

3. In the model, the accretion rate depends on the 
radial velocity and the surface density. Although radial 
velocity increases, but surface density reduces in the pres- 
ence of wind. As we analyzed the solutions of a and j3 
disks, the accretion rate reduces as the wind becomes 
stronger. Reduction to the accretion rate is not very 
sensitive to the value of s, but parameters I and Ao have 
a more significant effect to this reduction. 

4. Since the structure of the disk is modified in the 
presence of wind according to the solutions, we also stud- 
ied gravitational stability of the disk via Toomre param- 
eter. As the wind gets stronger, we see that Toomre 
parameters becomes larger which implies a more stable 
disk. However, dependence of Toomre parameter to the 
wind parameters are not at the same level. For example, 
Toomre parameter is not very sensitive to the exponent s. 
As it is shown in the figures, the gravitational instabili- 
ties in /3-disks are more pronounced than a-disks, even in 
the presence of wind. So, it might be anticipated that the 
/3-model can better describe the planet formation around 
new-born stars. In the case of proto-planetary disks j3- 
prescription yields the spectra that are considerably flat- 
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ter than those due to non-self-gravitating disks, which 
is in a better agreement with observations (Abbassi & 
Ghanbari 2008). 

5. As may be inferred from the figures, /3-disks typi- 
cally have an aspect ratio H/r smaller than a-disks, and 
thus fall into the thin disk regime to a greater degree. 
However, this ratio increases as wind get stronger, irre- 
spective of the viscosity prescription. 

6. We also found that thermodynamics of the disk 
has a vital role even in our simplified picture in which 
a poly tropic equation of state is used. With increasing 
7, there are reductions to radial velocity, surface density, 
accretion rate and mass loss efficiency, but the angular 
momentum per unit mass increases. Moreover, the disk 
is geometrically thicker and gravitationally more stable 
as 7 increases. 

The differences between a and (3 models of viscosity 
prescription were predicated by Duschl et al. (2000), and 
indeed is confirmed by our results. In a global overview, 
as in AGS06, we have shown that in the outer part of 
the disk, where the self-gravity has an influential func- 
tion, these models behave differently. They are somehow 
similarly affected by wind, though. In the real accretion 
disks, there are several important processes other than 
viscosity and wind, which are also expected to trans- 
port angular momentum outward. It is also immediately 
clear that the changes in the boundary conditions affect 
the structure of the solutions. Many questions remain 
about the wind itself. For example, how is it driven and 
where does it leave the disk? We, however, did not con- 
sider the driving mechanisms of the wind. Given these 
facts, the treatments in the paper are very simplified, but 
sufficiently general to describe many of the disk- wind sys- 
tems. All told, we believe that in order to obtain a better 
physical picture of such systems, more careful treatment 
is required and the analysis must be as deep as possible. 
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Fig. 3.— The profiles of the physical variables for a' = 0.1, t = 1.0, A/ in f a ii = 2.0 X 10 -6 M Q yr _1 and I = 0.0, 0.5, 1.0, 1.5 with A = 0.1 
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labeled by corresponding (. No- wind solution is shown by dashed curves. 









6 r- 







4 - 







2 - 


>° 





- 




-0 


2 - 


>° 


-0. 


4 - 








>^ 


-0 


6 - 




-0 


S - 




S 3 
s 

~"o 2 
g 

5 1 





0. 8 r 




0.7 - 




0.6 : 


- s 


0.5 - 
0.4 - 




i 0. 2 = 



10 



m 


05 p 





04 - 


" . 


03 - 


O 




V 


02 - 


- 


01 - 


3 n 


00 - 





0. 












• - : - " 



10 ■ 





10 3 10 2 10' 10° 















1. 5 












1.0" 












-0.0^= 















10" :! 10" 2 10" 1 10 



Fig. 4.— The profiles of the physical variables for /3' = 10 -3 , M in f a ii = 2.0 X 10 -6 M yr" 1 and I = 0.0, 0.5, 1.0, 1.5 with A = 0.1 and 
s = 0.7 at 7 = 1.1. Surface density and the rotational velocity for no-wind solution are represented by <jq and tio,^. Each curve is labeled 
by corresponding I. No-wind solution is shown by dashed curves. 
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